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ABSTRACT 

We confirm that the standard assumption of isothermal, shock-heated gas in cluster 
potentials is unable to reproduce the observed X-ray luminosity-temperature relation 
. of groups of galaxies. As an alternative, we construct a physically motivated model 

Ph ' for the adiabatic collapse of pre-heated gas into an isothermal potential that improves 

OO . upon the original work of Kaiser (1991). The luminosity and temperature of the gas 

fSJ ' is calculated, assuming an appropriate distribution of halo formation times and radia- 

tion due to both bremsstrahlung and recombination processes. This model successfully 
' reproduces the slope and dispersion of the luminosity-temperature relation of galaxy 

groups. We also present calculations of the temperature and luminosity functions for 
jy-^ ' galaxy groups under the prescription of this model. This model makes two strong pre- 

dictions for haloes with total masses M < 10 13 M Q , which are not yet testable with 
current data: (1) the gas mass fraction will increase in direct proportion to the halo 
mass; (2) the gas temperature will be larger than the virial temperature of the mass. 
The second effect is strong enough that group masses determined from gas tempera- 
tures will be overestimated by about an order of magnitude if it is assumed that the 
' gas temperature is the virial temperature. The entropy required to match observations 

can be obtained by heating the gas at the turnaround time, for example, to about 3 
x 10 6 K at z=l, which is too high to be generated by a normal rate of supernova ex- 
plosions. The isentropic model breaks down on the scale of low mass clusters, but this 
is an acceptable limitation, as we expect accretion shocks to contribute significantly 
to the entropy of the gas in such objects. 



Key words: galaxies: clusters: general - intergalactic medium 



1 INTRODUCTION 



The abundance of virialised clusters and groups of galaxies can provide a sensitive probe of the amplitude of density 
fluctuations on different scales and, hence, can be used to determine cosmological parameters. This requires, however, 
a theoretical framework for relating observable quantities such as the X-ray luminosity (L) and temperature (T) of 
the intra-cluster or intra-group medium (hereafter referred to as the ICM in both cases, for simplicity) to the total 
mass of the structure. The most commonly applied model assumes that the gas has been shock heated to the virial 



temperature of the isothermal mass. This results in a scaling relation L cx T (Kaiser 1991; Eke et al. 1996), which 



does not match the observed relation for clusters, (approximately L cx T 3 (Edge, Stewart 1991; Markevitch 1998)) 
or that of groups, which may be even steeper (Ponman et al. 1996). It is also true that the L cx T 2 scaling does not 



require isothermal gas; we develop such a model in Appendix |A|. 

It was originally shown by Kaiser (1991) that, if the gas was pre-heated (for example by supernova explosions) 
before falling into the dark matter potential, the observed cluster L-T relation can be more accurately reproduced. 
Assuming adiabatic collapse under these conditions, Kaiser's self-similar model predicts L cx T 35 . However, in this 
case, the radial extent of the gas depends strongly on the size of the potential; for low mass groups, the gas will 
extend well beyond the virial radius and, for high mass clusters, it will be concentrated in the centre. This may 
result in unreasonable gas mass fractions within the virial radius. 

Other authors have attempted variations on the self-similar model. Bower (1994) explored arbitrary scaling 
relations of the form L cx (1 + z) s M p , where s and p are free parameters, in an attempt to reduce the predicted 
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evolution of the X-ray luminosity function. A special case of this treatment is the model of Evrard & Henry (1991), 
which yields the relation L cx T 11 / 4 . Both of these models are based on the assumption that the core gas possesses a 
minimum entropy, the value of which is determined by the entropy of the initially pre-heated gas. Outside this core, 
it is assumed that the gas is distributed isothermally. However, the physical reality of a gas core has not yet been 
demonstrated, and the assumption of a largely isothermal profile may not be supported by observations ( Markevitch 



et al. 1998). Furthermore, as noted by Ponman et al. (1996), this relation, when extrapolated to lower mass scales, 



does not fit the X-ray data of galaxy groups. 

Recently, Cavaliere, Menci & Tozzi (1997; 1998), have constructed a model in which the observed L-T relation 
on both cluster and group scales can be reproduced by varying the gas density at the virial radius, according to the 
accretion-shock strength, as determined by the temperature difference between the infalling and virialised gases. In 
the present work, we explore the case in which pre-heated gas is assumed to collapse adiabatically, with the effect 
of shocks neglected. This is somewhat similar to the model of Cavaliere et al. on group scales, where shocks are 
expected to be weak, although we put physical constraints on the total gas mass within virialised haloes and do not 
require an isothermal distribution. 

The paper is organised as follows. In § our cosmological and structure formation parameters are defined. The 
specific models for the gas distributions are presented in § ^: in particular the isothermal model in § 3.1 and the 
isentropic model in § 3.2.2. The luminosity-temperature-mass relations and the temperature/luminosity functions 
are presented in § || The models are discussed in terms of their evolutionary predictions, energy requirements and 
failures in § |5|. The conclusions are summarised in § ^|. 



2 THEORETICAL FRAMEWORK 



2.1 Cosmological Parameters 

We will restrict our analysis in the present work to two simple cosmologies: a standard cold dark matter (CDM) 
model with fi Q = 1, h = 0.5, and an open model with f2 Q = 0.3, h = 0.75. The values of h are chosen so that the 
age of the universe is equal to about 13 Gyr in both models. Unless otherwise specified, the baryon fraction of the 
universe, fit,, will be given by the big bang nucleosynthesis value = 0.0125ft-" 2 (Copi, Schramm, Turner 1995). 

A simple function is fit to the z = variance of mass fluctuations, a(M,0), from the CDM power spectrum 
(Bardccn et al. 1986) for the two cosmological models: 



log(a(M,0)/a 8 ) = - 



0.0128 + 0.298 log 
0.017 + 0.248 log 



M 
Ms 
M _ 



f 0.0211(log 
0.0187(log^) 2 



M 



; = 1 , 
; n = 0.3 



(1) 



where the normalisation of the power spectrum is determined by the constant a% = a(Ms), and Ms = ^-p c (0)fl o (8h~ 1 Mpc) 
5.962 x lO 14 /i _1 r2 o M0 is the mass within an 8/i _1 Mpc sphere, with p c (0) representing the critical density at z = 0. 
These fits correspond to a spectral index n (P(k) cx k n ) ranging from n = —2 at the low mass end to n = —1 for 
the largest masses considered in this work. 

The linearly extrapolated value of <?(M, 0) at a redshift z is <r(M, z) = a(M,0 )g(z)/g(0), where the growth term 



g(z) is given, to a close approximation in the case of fl < 1 ( Viana, Liddlc 1996 ), by 



g(z) = (l + z)- 
where Sl(z) = o (l + z)/(l + Q Q z). 



n(z) 

2 



n(z) 



4/7\ 



Sic 



= i, 
< i, 



A = 



(2) 



2.2 Halo Formation Times 

For a halo of a given mass, the density profile of the ICM is related to the shape and the depth of the gravitational 
potential, with the latter depending on the epoch of halo formation; haloes of a given mass that form at an earlier 
epoch are denser, and more compact. This effect has often been neglected in the past, but has been treated recently 



by Kitayama & Suto (1996a, 1996b). Numerical simulations (Navarro, Frcnk, White 1995) suggest that the depth of 
the potential well, as traced by the circular velocity V c , remains relatively unchanged after 75 per cent of the cluster 
mass is in place, since the rest will be accreted in minor merger events that do not significantly disrupt the mass 
already in place. Thus motivated, we will define the redshift of halo formation as the redshift at which 75 per cent 
of the mass has been assembled; the X-ray temperature of the gas corresponding to a halo of mass M at the epoch 
of observation will then be determined by the circular velocity of that halo when its mass was equal to 0.75M. We 
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Figure 1. The probability distribution of halo formation times in the mod el o f Lacey & Cole (1993), as parameterised by lo, which is 
related to the halo mass and redshift of formation as described in the text (§ 2.2). The three different curves correspond to three different 
definitions of halo formation; the factor / represents the fraction of the final mass (at the observed epoch) that is virialised at the epoch 
of formation. 



will assume that both gas and dark matter components of merging haloes mix thoroughly and that the mixed gas 
immediately settles into hydrostatic equilibrium after a merger event. 

An analytic description of halo formation times, applicable to any bottom-up hierarchical model in which struc- 
ture formation is due to gravitational instability, was developed by Lacey & Cole (1993; 1994). For a halo of mass 
M, observed at redshift z, the probability that a fraction / of its mass was virialised at an earlier redshift Zf is 



M z f) 



2u(f- 1 - l)erfc(^L) - 



du) 
dzf 



■(f- 1 -2) — , (3) 



where the parameter u is given by: 



u(z f ) 



5(0) 



f) 



i) 



9(z) yj a 2 (fM(z), 0) - a 2 (M(z),0) 

1/2 



(4) 



Equation || is only strictly true for the special case a(M,0) — a s yj^J , but the dependence on the shape of 

the power spectrum is weak (Lacey, Cole 1993) and so is neglected for simplicity. The dependence on the parameter 
/ is shown in Figure [j] the most probable value of u> decreases from 0.6 to 0.3 as / increases from 0.65 to 0.85, 
but this has a negligible effect on the present calculations. The distribution of formation times, assuming f=0.75, 
is shown in Figure ^ for three representative mass scales at the present epoch, in both cosmological models. Lower 
mass haloes have a broader distribution of formation epochs, and haloes in the low density universe model form at 
earlier times, on average, than those in the critical density model. This analytic description of halo formation times 
has been shown to agree well with numerical simulations ( |Lacey, Cole 1994 ; Eke, Navarro, Frcnk 199§| ). 



2.3 Construction of the Temperature and Luminosity Functions 

The comoving abundance of clusters as a function of mass can be determined using the Press-Schechter (1974) 



formalism (also Bond et ah 1991), which has been supported by some numerical simulations (Efstathiou, Bond 
White 1992| ; |White, Efstathiou, Frenk 1993| ; [Lacey, Cole 1994| ; [Eke ct al. 1996| ) and is given by 



dN(> M) dM = 12 Pc (o)nj, 



dM 



exp 



na(M,z) 2 M ^\2a(M,z) 



da(M, z) 
dM 



dM, 



(5) 



where 5 C = 1.686fi(z) 0185 is the value of the linearly extrapolated critical overdensity for spherical gravitational 
collap se (corresponding to the top-hat fi lter) at redshift z, assuming no cosmological constant QNavarro, Frcnk, White 
1996 ) . There is accumulating evidence ( Bertschingcr, Jain 1993 ; Sovcrnato et al. 1998 ; gomervillc et al. 1998|) that 



the shape of the mass function on the scales of groups and clusters, as determined from higher resolution numerical 
simulations, is not the same as that predicted by the Press-Schechter formalism. The differences amount to a factor 
of 2-3 on both high and low mass scales; nonetheless the Press-Schechter prediction will be used here, and a study 
of the error introduced by this procedure will be deferred. 
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Figure 2. The probability distribution of halo formation times as a function of redshift in the model of Lacey & Cole (1993), for haloes 
at z = with three different masses, assuming f=0.75. The top panel shows the results for Qo = 1, while the bottom panel shows the 
results for Q = 0.3. Compared to high mass haloes, the low mass haloes have a higher most probable redshift of formation, as well as a 
broader distribution of formation times. 



The differential temperature function is obtained by integrating the mass function, Equation over all redshifts: 

dN(> T) _ dN(M(T, Zf )) dM(T, zj) dp{z f ) 

dT ~ J z dM dT dz f h [ ' 

where z max is the redshift at which Comptonization will efficiently cool virialised haloes. The value of z max = 7 
is adopted, but the results are insensitive to this choice for reasonable values (i.e., z max > 5). The cumulative 
temperature function N(> T) is obtained by numerically integrating Equation || over T. Thus, a relationship 
between mass and temperature is required, and this will depend on the thermal history of the ICM. An analagous 
procedure is followed to derive the luminosity function, given a relationship between luminosity and mass. 



3 MODELS 

A singular, truncated isothermal sphere will be assumed for the dark matter potential, p(r) — Pn(r / R v i T )~ 2 , where 
Pr is the density at the virial radius, i? v in and is equal to a third of the mean density within R vlr , p(R vlI ). This latter 
quantity is related to the critical density at redshift z by p(R v [ r ) = A c (z)p c (z), and A c = 78fl(z) + 80 + 300£!(z) /(l + 
I5ft(z)) is a fit, accurate to better than 2 per cent, to the results of the spherical collapse model as presented in Eke 
et al. (1996). It will be convenient to define a redshift evolution term, F(z) 2 = (1 + z) 2 (l + Cl z)A c (z)/ A c (0), so 
that 

p R = iA c (0) Pc (0)F(z) 2 . (7) 

For fl = 1, F(z) 2 = (1 + zf and A c = 178. 

The mass M of a cluster at its formation redshift Zf will be defined as the mass within i? V j r . Thus, the circular 
velocity V c , which is independent of r, is given by V 2 = GM/R vlr . The virial temperature is related to V c by 

kT vil = ^pm H V 2 (8) 

= 3.1( -| k.-Y. 

V1000 km s" 

From M — AttprR^. the following relations are obtained: 

M = V^Anpn)- 1 /^- 3 / 2 



A c (0)Y 1/2 ( V c 0) 



2.46 x wPh-^Fiz,)- 1 ^ ' e M G 



178 / V 1000 km s 



© 0000 RAS, MNRAS 000, 000-000 



Isentropic Gas in Galaxy Groups 5 



and 

kT vil = 0.5Gpm H {ATrp R ) 1/3 M 2/3 

'A c (0)\ 1/3 , , 2/3 f M \ 2/3 , (10) 



The volume emissivity of the X-ray emitting gas at radius r is given by 

e(r) 3_*L*IkM (n) 
2 /jm H t coo i 

where T g (r) is the gas temperature, and i coo i is its cooling time. For the latter, a fit ( Peacock 1996[ ) to the cooling 
function of Raymond, Cox & Smith (1976), constructed for an H and He plasma with Y = 0.25 (// = 0.59) and 
incorporating cooling due to both bremsstrahlung and recombination, is implemented. Recombination radiation is 
important for haloes with kT < 4 keV, and its inclusion allows us to consider the mass regime of galaxy groups. The 
expression for the cooling time is 

<cool(r) = ^frMWl (12) 

where C\ = 3.88 x 10 11 s K" 1 / 2 cm" 3 , C% = 5 x 10 7 K and f m is a metallicity dependent constant that is 1.0 for 
solar metallicity (adopted here), and 0.03 for zero metallicity. We only consider haloes in which not all of the gas 
within i? v i r has had time to cool since the halo formed. The luminosity is given by integrating the volume emissivity 
out to the virial radius: 

"-Rvir 

Airr 2 e{r)dr 

R ( 13 ) 

" C^f I ' 'WWII + ft/.T. W -'l*. 

The final density profile of gas that is accreted by a dark matter halo depends on the density profile of the dark 
matter itself and also on the difference between the entropy of the virialised and infalling gas components. If the 
entropy of the infalling gas is low, it will be heated by accretion-shocks to approximately the virial temperature of 
the halo, T v ir; thus, for our adopted mass profile, the gas distribution will be isothermal. Simulations show that, 



indeed, cluster temperature profiles are isothermal to within a factor of 2-3 (Navarro, Frenk, White 199£; Eke 



Navaro, Frenk 1998). Alternatively, if the gas is pre-heated so that it has a large initial entropy, it will be accreted 
approximately adiabatically, and the gas profile will be isentropic. The purpose of this paper is to consider the 
properties of cluster gas under these two extreme scenarios. For the majority of clusters, the true physical process 
is probably intermediate between these models, with the gas becoming more isothermal and less isentropic with 
increasing halo mass. 



3.1 Isothermal Model 



In this model, we make the common assumption (e.g., Eke et al. 1996) that the gas is distributed isothermally, 
with a temperature equal to the virial temperature of the halo. We include a distribution of formation redshifts and 
recombination radiation, which are often neglected in studies of this nature. The formation redshift Zf of a halo 



of observed mass M has been defined (§ 2.2) to be the redshift at which the mass was equal to fM. Since the gas 



temperature is assumed not to change between zf and the epoch of observation, the observed temperature and mass 
are related, from Equation by 

1/3 / f \ 2/3 / , , v 2/3 



kT n = 4.5817 



2/3 



A c (0) 



/ 



2/3 



M 
Ms 



keV. 



(14) 



178 ) V - 75 , 

This gives lower temperatures than the relation of Eke et al. (1996), due to the factor of f; it is only slightly lower 
than the normalisation of Pen (1998) which was determined from numerical simulations. 

If the gas is dissipationless, its density profile will match that of the dark matter, i.e., p g (r) = p gi R(r/i? v ir)~ 2 , 
and Pg,R,/ Pr = Slb/fio- To avoid the singularity at r — when integrating over the assumed isothermal profile, an 
arbitrary core radius of r c = f c R v ir is adopted with f c — 0.1, such that p g (r < r c ) = p g (r c ). 

Evaluating Equation [Tj] yields the following luminosity-temperature relation: 



L = 5.83 x I0 42 h 



f2j,/fi c 
0.05 



A c (0) 
178 



[f) \keVj \ 12.33 



1 + 4.31/„ 



^keVy 



ergs 



(15) 

The luminosity can be expressed as the sum of two components ibrem and L rcc , corresponding to the bremsstrahlung 
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and recombination radiation, respectively. Substituting T g from Equation |l4|, Lbrem oc M 4 / 3 i 7 '(zy) 7 / 3 and L r 
M 2 / 3 F(z f ) 5 / 3 . For n = 1, this corresponds to i b rcm oc A/ 4 / 3 (l + z f ) 3 - 5 and L ICC oc M 2 / 3 (l + z f ) 2 - 5 . 



3.2 Pre-Heated Gas 

3.2.1 Principles and Definitions 

The equation of state of a polytropic gas is P = K pJ, where P is the gas pressure, p g is its density, and 7 and K 
are constants. K Q is related to the specific entropy of the gas, S, and the specific heat capacity at constant volume, 
Cv, by ([Sears, Salinger 1975 ) 

K °= 2n(X H )^ eMS/Cv - 5m - (16) 



In the isothermal model (§ B.l), it was assumed that the gas was shock heated to the virial temperature upon 
accretion onto the dark matter potential; in this case, 7 — 1. At the other extreme, which is the focus of this section, 
gas with high entropy will be accreted adiabatically, and 7 = 5/3. Evaluation of K Q will be discussed further in 
5 13.24 



The gas density profile is obtained by solving the equation of hydrostatic equilibrium: 

dP _ GM(r) 
dr 



-Pgir)- 



(17) 



In principle, M is the mass of both the dark matter and the gas, but we will neglect the contribution of the gas, the 
mass of which is only a few percent of the total (e.g., White et al. 1993). The resulting gas density profile within the 
virial radius is 



Pair) = p g , R [l + 3A- 1 ln(i? vir /r)]— ; 
A= 9c^(ii vir )_37^R 1 



(18) 

2 V 2 (7"1)K 2 ' (19) 

where c s (i? v i r ) = \f^K p^ K is the adiabatic sound speed at the virial radius. This density gradient is significantly 
shallower than that of isothermal gas and, as a result, the X-ray surface brightness profile is quite flat near the centre, 
producing an X-ray core without requiring a core structure in the underlying dark matter or gas mass distributions. 
Assuming 7 = 5/3, we derive the ratio of gas mass to dark matter mass within the virial radius from Equation 18: 



= 0.33^A- 3 / Vr(2.5, A), 



M p R 
~ t t a ~ 1 dt is the incomplete gamma function. 



where T(a,x) = J c e 

For an ideal, polytropic gas, the temperature is related to the density by 



Pa 



( kT g 



\p:mHKc 



so 



where 



T g (r) - T g>R [l + 3 A- 1 ln(i? vir /r)] 



kT g , R = 7 l pm H V 2 A = |- — -AkT vir 
37 3 7 



(20) 



(21) 



(22) 



(23) 



is the gas temperature at the virial radius. Observationally, the mean cluster temperature measured is an emission 
weighted temperature, given by 

/ e{r)T g (r)r 2 dr 



T — 



The relationship between T cm and T 3j r is 
kT - 



fcT 9jfl r(5.5,A) 



A r(4.5,A) 



J e(r)r 2 dr 
4.3L4/™ T(4.5, vl) 



4.3LA/ m r(3.5,A) 



kT g , R T(4.5,A) 



(24) 



(25) 



kT g>R r(5.M)_ 

The bolomctric luminosity is given by substituting the density and temperature profiles from Equations 18 and 
into Equation [l3| The result is (assuming 7 = 5/3): 



L 



6.3 x IQ-^k 1 /^} 



[C 1 p 1 m H G^ c {0)pM] lb K 3 o F{zf 



e A r(4.5,^) 



7.5fcC 2 / m r(3.5,^) 
fim H V 2 r(4.5,A) 



(26) 
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3.2.2 Relating K a to the Gas Temperature 

Here we will assume 7 = 5/3. The entropy constant K 0l denned in Equation |l6|, is then related to the temperature 
and density of the gas by 

kT 

K ° = hi* ( 2? ) 

p,m H p g 

-2/3 

/ f . \ I I J . . \ 

= 3.70^34 



( T g \ ( Pg 



VlO fi K; \Q bPc (0), 

where K34 = 10 34 ergs g -5 / 3 cm 2 . 

The entropy of pre-collapse gas can be increased by heating it to a temperature which will depend on the gas 
density at the time of heating. We will consider heating the gas at the turnaround time, when the overdense region 
decouples from the Hubble expansion and begins to collapse; if the heating takes place instead at some later time 
during the collapse of the halo, the gas density will be higher and the amount of heat ing necessary will be greater. 
At the turnaround time, the radius of the halo is equal to twice the virial radius e.g. (Lacey, Cole 1993) and, thus, 
if we assume the perturbation to be homogeneous, we have 

Pta = gpORvir) = ^A c (Q) Pc (Q)F(z f ) 2 . (28) 
Therefore, we can relate the entropy to the gas temperature at turnaround, T g ta , by 

K = 314.8^34 [A c (0)F(z } Y}- 2/3 - (29) 

In the special case of tt = 1, Equation ^9] becomes 

K = 0.47^34 (1 + zfY 2 , n = 1. (30) 

The temperature required to obtain a given entropy decreases for lower values of f2 , due to the o dependence of 
A,, 



3.2.3 The Isentropic Model 

Assumptions need to be made about the behaviour of p Sl R and K a to evaluate the parameter A in the equations 



of § 3.2.1. A simple but physically unmotivated assumption is the one we made in the isothermal model, that 
Pg.R./ PR = ^b/^o- We discuss interesting results of two such models in Appendix^, but neither produces an entirely 
satisfying result. Furthermore, constraining p gi R at r — 00 is not a satisfactory constraint because gas cannot be 
assumed to be in hydrostatic equilibrium beyond a distance approximately equivalent to c s /H , where H a is the 
Hubble constant; this distance is equal to a cluster virial radius within a factor of order unity. Beyond this radius, 
gas will be infalling and its distribution is not represented by our model. 

Instead, we will assume that the entropy constant K is independent of halo mass, and that the total mass of 
gas in the cluster is determined by the amount that could have been accreted at the adiabatic Bondi accretion rate 



([Bondi 1952ft , assuming that the accreted gas has already been pre-heated, and accretes isentropically: 

M = 4^AG 2 M 2 ( 7 ^or 3/2 pf (5/3 " 7) , (31) 
w l.86n\G 2 M 2 Ko 3/2 , 7 = 5/3 

where A = 0.25 is the dimensionless accretion rate and p g and T g are the density and temperature of the gas being 
accreted. We assume 7 = 5/3 for the remainder of the derivation. If the simplifying assumption is made that the 
accretion rate has remained roughly constant, the amount of gas mass accreted by the observed redshift z, M g {z), 
can be determined from 

^# = 0.040^ ^ , (32) 

M 10 14 M Q \K 34 J 10 9 years' V ' 

where t(z) is the Hubble time at redshift z. Thus, we predict a strong mass dependence of the gas mass fraction, 
Mg/M ex M, for the lowest mass haloes (see below for the mass limit). There is currently insufficient data available 
to accurately test this prediction, but there are indications that the gas mass fraction increases significantly from 



elliptical galaxies to groups to clusters (David, Jones, Forman 1995) 



The model gas mass fraction, from Equation |32|, will exceed Ob/Op for clusters with circular velocities greater 
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than the critical value V^ r it : 

. — (^>f (^f (±f (j^y lf °^ ta .-. ,33, 

For Q = 1, this corresponds to a halo mass of M crit = 3.78 x l0 14 hil b (K o /K 34 )^ 2 (l + z f fl 2 M & . Above this 
mass limit, the Bondi accretion rate will be ignored, and the condition M g /M = f2ft/f2 will be imposed. This is in 
reasonable agreement with observations (e.g. Edge & Stewart 1991; Evrard 1998) which show that gas mass fractions 
are roughly constant among clusters, though there is a considerable amount of scatter. The critical halo mass should 
be compared with the Jeans mass at the turnaround time, which is defined by 

47T 

Mj = —tt oPc {z ta )\,{z ta f (34) 

where the Jeans length is Aj = t(z)c s . Using the relations tf = 2t ta and 4(1 + zj) 3 = (1 + z ta ) 3 , (the latter valid 
only for Q = 1), the Jeans mass of a halo at turnaround can be expressed in terms of its formation redshift 
as Mj = 2.4 x lQ 1A hVt b {K /K 34 f/ 2 (\ + z f f/ 2 M Q . Thus, we obtain M clit /Mj = 1.6 which is independent of of 
Zf, K and h. Since the Jeans mass, which determines the scale at which pressure effects are important relative to 
gravitational collapse, is not much smaller than Mb, our assumption of adiabatic Bondi accretion is not unreasonable. 

For low mass haloes, with circular velocities less than the critical one given in Equation [33], the constant p g .R is 
defined parametrically (through A) by 

^,A ) = e ^^) 1,2 F M J^- a . (35) 

Thus, the parameter A is independent of V c , but does depend on the redshift of formation in a non-trivial manner. 
From Equation 19, this requires p g .R oc V 3 . Since A is independent of V c then, from Equations ^|, 23 and |25| 
the isothermal M-T relation T cm oc M 2 / 3, is again obtained, but the constant of proportionality changes and the 
dependence on formation redshift is more complex. Allowing for the fact that, on average, higher mass haloes form at 
lower redshifts than low mass haloes, the M-T relation will be very slightly modified. The luminosity dependence on 
mass, however, being more sensitive to the gas distribution, differs significantly from the isothermal case. Considering 
bremsstrahlung radiation only, for example, L oc V^°F(z)~ 3 oc T e 5 m F(z) -3 . Thus, the important result is obtained 
that, by assuming adiabatic collapse of gas and a gas mass fraction which is proportional to the total mass, L oc T 5 m . 

For high mass haloes, with circular velocities greater than the critical value, the situation is considerably more 
complex. To solve /0 gi R, the parameter A must be determined from the following equation: 

M.5 (Mb 



e A T(2.5, A) = (7.5K Q )^ (-£ j [A c (0) p c (0)]^ 3 F(z) 2 . (36) 

Thus, A (and, hence, p g ,R) depend on both the mass and formation redshift in a complex manner. Furthermore, 
Equation Bra is not soluble for all V c . The left hand side of the equation approaches 1.33 as A approaches zero, which 
means that solutions can only be found for haloes with circular velocities that satisfy 

The model breaks down because the gas becomes increasingly concentrated toward the cluster centre as mass increases 
and, eventually, the density at the virial radius becomes zero. Since shock heating is expected to dominate the gas 
entropy for large masses, this is an acceptable limitation. Alternatively, a solution exists for any mass if the entropy 
constant K is made sufficiently high. 



The bolomctric luminosity, L, is then determined from Equation |26| after solving for A from Equation |35| or 36 
Similarly, the emission weighted temperature is computed from equations |23| and ^5|. Note that L and T do not scale 
simply with h, due to the complex h dependence of A. 

The restriction on circular velocity presented in Equation ^ can be expressed in terms of temperature and 
luminosity limits as: 

and 

. 11 /o / K. "n \ „ , s 11 /q /fl/i/i^n \ { f ( 1 ^ 

(39) 



L crit < 1.61 xl0 44 /i n / 3 — F(z) 11 / 3 

1 K 34 J 1 ' V 0.05 J V 178 



K 34 J \ 0.05 178 



ergs s . 



© 0000 RAS, MNRAS 000, 000-000 



Isentropic Gas in Galaxy Groups 9 




.o 



g T w [ergs s l ] 



Figure 3. The luminosity-temperature relation, with data from David et al. (1993,1995 stars), Markevitch (1998, circles), Ponman et 
al. (1996, open squares, resolved sources only) and Mulchaey & Zabludoff (1998, open triangles, Raymond— Smith model determinations 
with half-solar metallicities). The error bars are omitted from the Markevitch data, as they are smaller than the data symbols. The 
isothermal model and isentropic model are plotted; the hatched region reflects the range of formation redshifts for clusters of a given 
mass. Only the isentropic model matches the slope and dispersion of the group data well; the normalisation of this model depends on 
the initial entropy, which we determine by fitting to the group data. The straight, solid line represents the model of Evrard & Henry 
(1991). Although this matches the cluster data well, it underpredicts the temperature of low luminosity groups. The models and data 
do not scale with h in an identical manner, as discussed in the text. The difference between the top panel (Q a = 1, h = 0.5) and bottom 
panel (H = 0.3, h = 0.75) is due primarily to the different values of h used, and not Q - 



4 RESULTS 

4.1 The Luminosity— Temperature Relation 

The L-T relation]^ for both models and cosmologies is shown in Figure^; temperature is plotted on the vertical axis, 
somewhat contrary to convention, since this quantity has the dominant observational uncertainty. For each model, 
three lines are drawn, representing the range of formation redshifts between Zf = and the la upper limit (i.e., where 
the probability that the cluster formed at a lower redshift is equal to 68.3 per cent). The central line represents the 
calculation for the most probable formation redshift. The data plotted are those of David et al. (1993;1995, stars), 
Markevitch (1998, circles), Ponman et al. (1996, open squares )[J] and Mulchaey & Zabludoff (1998, open triangles)^. 
All data points are unique; where there is overlap we choose the Ponman et al. group data over the Mulchaey & 
Zabludoff data, and the Markevitch cluster data over that of David et al. 

It is crucial to note that the data and models do not scale with h in the same way; for this reason, h has not 



* The emission weighted temperature of the model is always considered when comparing with observations. 
T Only fully resolved observations are considered (i.e., with a quality index of 1). 

t We use the temperatures determined using the Raymond— Smith model with the metallicity fixed at half solar for all groups except 
NGC5846, for which this temperature is unconstrained. In this case, we adopt the low metallicity determination. 
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been left as a free parameter. The observed luminosities scale as L oc h~ 2 , while the isothermal model0 scales as 
L oc h£l%, from Equation [l5|. Thus, if fib is fixed, L oc h, while if the big bang nucleosynthesis value is adopted, as we 
do here, L oc h~ 3 . In either case, the ratio of model to observed luminosity is h dependent. This accounts for most 
of the difference between the matching of the model and the data in the f2 = 1 cosmology (top panel, h = 0.5) and 
the Cl = 0.3 cosmology (bottom panel, h = 0.75). 

The slope of the isothermal model is too steep to match the cluster data, even for haloes with temperatures greater 
than a few keV. Due to the effects of recombination radiation, it steepens even further for kT < 4 keV, approaching 
L oc T and, thus, completely fails to match the low-luminosity, group data. Furthermore, for luminosities less than 
about 10 43 ergs s _1 the gas will have had time to cool efficiently since the formation of the halo and will no longer 
be emitting X-ray radiation. 

For the isentropic model, the entropy constant, K , has been determined by matching the L-T relation to the 
available group data of Ponman et al. (1996) and Mulchaey & Zabludoff (1998). This constraint requires a value of 
K Q = (0.37 ± 0.1)K 3i which, from Equation |(| can be obtained by heating the gas at the turnaround time to a 
temperature of about T = 7.9 x 10 5 (1 + z/) 2 K, for Q = 1. For Zf w 1, the required gas temperature is 3.1 x 10 6 
K; this is reduced to T = 1.8 x 10 6 K for fl a = 0.3. 

The slope of this model for the lowest masses (L oc T 5 ) is certainly consistent with most of the group data. 
Ponman et al. (1996) claim the slope of their data is best fit by L oc T 8 - 2 , but there is a great deal of scatter and there 
are significant outliers. Including the data of Mulchaey & Zabludoff (1998) we obtain a best fit relation of L oc T 3 5 , 
but the chi-squared value indicates a poor fit, so the errors on the slope are not well determined. Clearly, more data 
with smaller errors are required to constrain the slope in this region. It must also be noted that the Ponman et al. 
(1996) sample is a collection of Hickson compact groups, which are atypical groups. It is, in fact, remarkable that 
they lie in the same area of the L-T plane as the Mulchaey & Zabludoff (1998) group data. 

The change in model slope at L ss 10 42 ergs s _1 is due to the change in gas mass constraint, from the amount 
of mass accreted at the adiabatic Bondi rate to a fixed value (at V c = V cl -it' 7 see § [3.2.3 ). As the halo mass is 



increased beyond this value, the central gas density increases rapidly, and the gas density at the virial radius drops. 
Since the density and temperature gradients are initially fairly shallow, the total emission weighted temperature is 
dominated by the greater volume of gas at large radii; the lower temperature of this gas approximately compensates 
for the central increase in temperature, and the total emission- weighted temperature remains roughly constant. 
The luminosity, on the other hand, is more strongly dependent on density and, thus, continues to increase with 
V c (although not as strongly as for haloes with V c < V cr it)- In this mass regime, the model does not fare as well, 
generally underpredicting the observed temperatures; this is not unexpected, as accretion shocks will begin to become 
important for haloes of this size. One possible model for the structure in these haloes is that the pre-heated gas 
dominates in the central regions, while accretion-heated gas dominates elsewhere. A simple model of this type 
is described by Evrard & Henry (1991), and is represented by logfcT = 4/11 \og(h 2 L) — 15.42 when normalised 
to the cluster data. This model is shown here as the straight, solid line, which our isentropic model approaches 
asymptotically. The Evrard & Henry model matches the slope of the cluster data quite well, but underpredicts the 
temperature of most groups and the lowest mass clusters. 



It has been shown (Fabian et al. 1994; Arnaud, Evrard; Markcvitch 1998) that the dispersion in the L-T relation 



can be reduced to a very small value (an r.m.s. dispersion of about 0.11 in log L at a given T) by excluding cooling flow 
regions from the observed X-ray data. The predicted (lcr) range in log L for a given T for high mass clusters under 
the assumptions of the isothermal model, due to the distribution of halo formation times, is 0.07-0.14, consistent 
with this data (though, of course, the slope is wrong). The isentropic model predicts a larger range of about 0.3 for 
low luminosity clusters and groups, and this range encompasses most of the group data shown in Figure ||. This 
suggests that the observed dispersion in the L-T relation may be primarily due to the distribution of halo formation 
times. 

4.2 The Mass- Temperature Relation 

The mass derived from the X-ray temperature of group-scale structures is strongly dependent on the thermal history 
of the gas, as shown in Figure ||, where the mass-temperature relation is plotted for the isothermal and isentropic 
models; as in Figure ^, the hatched region represents the dispersion in the relation due to the distribution of halo 
formation times. At low masses, M < 10 13 M Q , the gas mass corresponding to a given temperature is about an order 
of magnitude smaller in the isentropic model than in the isothermal model: the mass of haloes with kT < 0.4 keV will 
be greatly overestimated by assuming the gas has been accretion-shocked to the virial temperature. Furthermore, 



§ Recall, from § 3.2.3, that the values of L and T in the isentropic model do not scale with ft in a simple manner. 



© 0000 RAS, MNRAS 000, 000-000 



Isentropic Gas in Galaxy Groups 11 



1 

: n — 1 

_ Isentropic ^ 


1 1 1 1 — : 


. Isentropic ^ 


Jp7 Isothermal 



10 12 14 16 



log M/M 

Figure 4. The mass-temperature relations for f2 = 1 (top) and f2 = 0.3 (bottom). The two models shown in each panel are the 
isothermal and isentropic models discussed in the text. The hatched region reflects the range of formation redshifts for clusters of a given 
mass. On the scale of groups of galaxies, both models predict M oc T 15 , but the mass corresponding to a given temperature is an order 
of magnitude smaller in the isentropic model, compared with the isothermal model. 



the dispersion of the M-T relation is significantly larger in the isentropic model; thus it will be difficult to determine 
the mass of individual, low temperature structures to better than a factor of ten or so from the X-ray observations 
alone. In the isothermal case, for high temperature clusters, the relation is much tighter, and the mass corresponding 
to a given temperature can be determined to within a factor of two. 

The most readily observable, independent parameter related to the halo mass is the line of sight velocity dis- 
persion, ay. Assuming isotropic orbits, we use the relation V c = \/2a r to compare our model predictions with 
observations. We show, in Figure |], the relation between ay and gas temperature predicted by the isothermal and 
isentropic models for £1 Q = 0.3. The group data plotted are the same as in Figure ||; velocity dispersions for the 
Mulchaey & Zabludoff (1998) are taken from Zabludoff & Mulchaey (1998). For the cluster data, velocity dispersions 
(not available for all of the data) are taken from Fadda et al. (1996) and various sources within the compilation of 
White et al. (1997). In the rare cases where uncertainties are not quoted, we assume 20% errors. 

The lowest temperature group data generally lie well below the isothermal model prediction; however, there is a 
great deal of scatter, and insufficient data below kT = 0.4 keV to provide a convincing test of our isentropic model 
prediction. Furthermore, even the massive cluster data tend to lie toward higher temperatures at a given ay, which 
may indicate that the galaxy orbits are not exactly isotropic. The most discrepant point is the cluster AWM4, which 
seems to have an unusually low velocity dispersion of 119 km/s (Fadda et al. 1996) for its temperature of 3.7 keV. 
Nonetheless, the data suggest a strong drop in ay, relative to the isothermal model, for haloes with temperatures 
below about 0.8 keV, in accord with our isentropic model. The notable exception to this trend is HCG 15, which 
has kT=0.44 keV and a r =457 km s" 1 (Ponman et al. (1996). 

The difference between the isothermal and isentropic models in Figures ^ and |B| is somewhat sensitive to the 
mass accretion rate A we have assumed in Equation 31. A lower accretion rate (i.e. A < 0.25) will preserve the slope 
of the M-T relation, while altering the normalisation, as the total gas mass within a halo of a given mass will be 
lower than we have assumed here. In this case, the temperature will drop. Our choice of A = 0.25 is the critical 
value determined by Bondi (1951) for the case of spherical, adiabatic accretion by a point source; we are currently 
exploring a more realistic accretion model, and indications are that the proper A may be somewhat lower than we 
have used here. 



4.3 The Temperature Function 

The cumulative temperature functions for the two models are shown in Figure ra. The Cl = 1 cosmology is shown 
in the top panel, and the Q = 0.3 results in the bottom panel. The z — model is shown as the solid line, and 
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Figure 5. We show the relation between line-of-sight velocity dispersion and gas temperature for the isentropic and isothermal models, 
with Q = 0.3. The width of the hatched region reflects the range of formation redshifts; the model with the wider band is the isentropic 
model. The data sources for the temperatures arc the same as for Figure B sources of the velocity dispersions are described in the text. 



Isentropic 
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Figure 6. The cumulative temperature function for n = 1 (top) and Q a = 0.3 (bottom). The isothermal and isentropic models are 
shown at three redshifts, z=0 (solid line), z=0.3 z=0.5 (dashed lines). The z=0 models are normalised to the data (Henry & Arnaud 
1991, with errors from Pen 1998) at kT=3.4 keV. For the isothermal case, evolution at the high temperature end is strongly negative 
for f2 = 1 and only weak ly so for Qo = 0.3. The isentropic model predicts some positive evolution for objects with kT < 2 keV but, as 
noted in the text (§ p.2.3|), it breaks down for larger mass haloes. 



the dashed lines represent the results at redshifts of 0.3 and 0.5. The cluster data from Henry & Arnaud (1991), as 
analysed by Eke et al. (1996) and with la errors taken from Pen (1998), is plotted as the step function in this figure. 
This data consists of the 25 X-ray brightest clusters in the HEAO-1 A2 survey, as detected in the 2-10 keV band. 

The r.m.s. mass fluctuation in spheres of size 8/i -1 Mpc, as, is determined by requiring the cumulative temper- 
ature function at z=0 to match the observations at kT — 3.4 keV, where the errors are smallest. For the isothermal 
model, erg = 0.802 for f2 = 0.3, and as = 0.492 for il Q = 1. These values are lower than those found by Eke et al. 
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Figure 7. The differential luminosity functions for Q = 1 (top) and Q = 0.3 (bottom). The two models shown in each panel are 
the isothermal and isentropic models discussed in the text, at z=0 (solid line), z=0.3 and z=0.5 (dashed lines). Very little evolution 
is exhibited in either model or cosmology. The flattening of the isothermal luminosity function at low L results from efficient cooling 
of the gas in low mass haloes. The data are from Rosati et al. (f998, open squares) and Henry & Arnaud (1991, filled squares). The 
normalisation of the models relative to the data is sensitive primarily to the value of h, since the two do not scale the same way. 

(1996), primarily as a result of including the distribution of formation times. Neglecting this (i.e., assuming each 
cluster formed at the redshift at which it is observed) overestimates as by 10-20 per cent, depending on cosmology 
(a larger effect for fl = 0.3). The effects of using a slightly different point of normalisation and mass-temperature 
relation are much weaker, affecting the results by about 1 per cent. 

The isentropic model is shown in Figure ^ for the same three redshifts. For this model, as cannot be calculated 
in the same way, as the model breaks down for temperatures greater than 1-2 keV, so the values determined above 
for the isothermal model are used. In this case, evolution is mildly positive; i.e., there are more objects of a given 
temperature at z=0.5 than at z=0, by a small amount that is roughly independent of cosmology. Relative to the 
isothermal temperature function, the isentropic model predicts about ten times more objects with temperatures less 
than 1 keV. This is also evident from Figure |^: the masses of haloes with a given X-ray temperature are much smaller 
in the isentropic model, and less massive objects are more numerous. 

4.4 The Luminosity Function 

Although temperature data do not yet exist for a complete sample in the low mass regime, a differential luminosity 
function is available from the ROSAT Deep Cluster Survey (RDCS, Rosati et al. 1998) to luminosities as low as 
2 x 10 42 ergs s _1 . We correct these 0.5-2 keV band luminosities to bolometric luminosities by dividing them by the 
factor exp(— 0.5/fcT) — exp(— 2/kT), where the temperature is estimated from the Evrard & Henry (1991) relation 
shown in Figure [| The lowest redshift data from the RDCS (0.045 < z < 0.25) are shown as the open squares in 
Figure [?]; the solid squares plotted on this figure are from Henry & Arnaud (1991). The isothermal and isentropic 
models are shown for both cosmologies. For each model, the solid line is the z = model, while the dashed lines 
represent the results at z = 0.3 and z = 0.5; very little evolution is exhibited in either cosmology by either model. 
This is consistent with the RDCS survey results, which show no evolution of the luminosity function out to z=0.8. 

The slope of the luminosity function corresponding to the isothermal model is too steep to match the data in 
the high luminosity regime, especially for the o = 1, h = 0.5 cosmology. As discussed earlier, the normalisation 
difference between the models and the data is mostly sensitive to h. The flattening of this slope at low luminosities 
is due to the fact that we do not include haloes in which the gas has had sufficient time to cool since the formation 
of the halo; efficient cooling in low mass objects restricts luminosities to those greater than about 10 43 ergs s _1 . 

The isentropic model, on the other hand, reproduces the flat slope (d\og(dN / dL) / d\og(L) = —1.88) of the 
RDCS data well, though the normalisation is too high. However, the major success of the L-T relation predicted 
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Figure 8. The evolution of the luminosity— temperature relation is shown for Q = 1 (top) and Q a = 0.3 for both isothermal and 
isentropic models discussed in the text. The solid line represents the relation at z=0, and the dashed line the relation at z=0.5. The 
(mild) evolution is in opposite directions in the two models. 



by this model occurs at luminosities lower than those probed by the RDCS (L < 10 42 ergs s _1 ) where the gas mass 
fraction is proportional to the total mass (§ 3.2.3). We predict that the slope of the luminosity function in this mass 
range will flatten to d\og(dN/dL)/d\og(L) = —1.33. 



5 DISCUSSION 



Recently, Henry (1997) published a detection of mild negative evolution of the cluster temperature function out 
to z=0.3, particularly at the hot end, kT > 5 keV. This is marginally inconsistent with the isothermal model in 
both cosmologies in Figure |[ as the fl = 0.3 temperature function shows too little evolution, while the evolution 
predicted by the fi = 1 model is too strong. For this reason, Henry (1997) claims f2 = 0.5 fits the data best, 
but there are few points and the error bars are large. The isentropic model predicts very little evolution in both 
cosmological models, but it is only valid at low temperatures for which there is currently no data. Since the redshift 
evolution in this case has such a weak cosmological dependence, estimations of cosmological parameters from X-ray 
abundances will best be made from the number evolution of the most massive clusters, and not from the abundance 
of lower mass clusters and groups. It is still important, however, to firmly establish the thermal history of the ICM 
before observed evolution in X-ray luminosity or temperature can be used to determine cosmological parameters. 

The direction of evolution of the L-T relation, shown in Figure ||, depends on the gas physics involved. In the 
isentropic model, clusters of a given temperature are significantly less luminous at z=0.5 than they are at z=0, 
whereas, in the isothermal case, they are somewhat more luminous. The work of Henry (1997) shows little or no 



evolution in the L-T relation for 3 < kT < 10 keV, in confirmation of earlier work (Henry et al. 1994; Donahue 1996 



Mushotzky, Scharf 1997). The prediction of moderate evolution in the isothermal model was part of the impetus 
for considering models with pre-heated gas ( |Kaiscr 1991 ; Evrard, Henry 1991); unfortunately, again, our isentropic 



model is not valid at these temperatures. Since errors in temperature tend to dominate observations, it is unlikely 
that the predicted evolution of the L-T relation on group scales will be easily observed. 

Many observations have suggested that the contribution of X-ray gas to the total mass in galaxy clusters is 



10-20 per cent, significantly larger than th e Big Bang nucleosynthesis value assumed here (White et al. 1993; White 



Fabian 1995 ; White, Jones, Forman 1997 ). Increasing the baryon density of the universe allows the application of 
our isentropic model to more massive haloes, as the limits in Equations |33| and |37] are increased. We exp lore, in 
Figure^, the luminosity-temperature relation assuming fib/fio = 0.1. (Below our critical mass, Equation the gas 
mass fraction is still given by Equation [3^ and, hence, less than f2b/f2 = 0.1, in qualitative accord with observations 
(David, Jones, Forman 1995)). This requires a slight modification of the entropy constant to K a = O.35K34 for 



© 0000 RAS, MNRAS 000, 000-000 



Isentropic Gas in Galaxy Groups 15 




o = 1 and K = O.39A34 for Q Q = 0.3. The improvement of the fit to the data is quite noticeable in the o = 1 
case, and is due primarily to the increase of the upper mass limit at which the gas fraction reaches Qb/fl . 

The isentropic model requires that the infalling gas have an initial entropy given by K Q — (0.37 ± 0. 1)^34. 
This is consistent, with about la confidence, with the observed value recently published by Ponman, Cannon, 
Navarro (1998), K — (0.20 ± 0.1)A"34. To generate this entropy requires some mechanism for heating the gas (at 
the turnaround time) from the background temperature of 10 4 K to about 3 x 10 6 K (see § 4.1); the energy required 
to heat a total mass of gas Sl^M from « 10 4 K to this temperature is 1.2 x 10 61 (£li,M/10 13 M@) ergs. If the amount 
of energy of a single SN event is about 10 51 ergs, 10 per cent of this is converted into thermal energy of the hot gas 
(Babul, Rces 1992), and the mass in stars is about 10 percent of the mass in gas, (White et al. 1993), this requires a 
SN rate of 1.2 x 10 9 events in the lifetime of a galaxy with stellar mass 10 10 M Q . Over the likely timescales of galaxy 
formation, this is probably an unreasonably high rate. 

Both models explored here fail to reproduce the temperature function, luminosity function and L-T relation for 
the majority of massive clusters. This mass range is successfully described in some aspects by the model of Evrard 
& Henry (1991), in which only the core gas retains its initial entropy, and the possibility of mergers stirring up 
and mixing the gas is ignored. An alternative model which may be at least equally successful is one in which the 
gas is more thoroughly mixed, and behaves like a polytrope with 7 varying from 5/3 on group scales to 1.1-1.2 on 
cluster scales, and approaching unity as the systems grow increasingly massive. Such a model has been postulated 
as an explanation of numerical simulation results (Lewis et al., in preparation), observed temperature gradients 
( Markevitch et al. 199S) and the discrepan cy between X-ray determined masses and those obtained from lensing 
observations ( Miralda-Escude, Babul 1995 ). We are currently exploring the predictions of such a model under the 
present formalism. 



6 SUMMARY AND CONCLUSIONS 

Two extreme models of hot intra-cluster and intra-group gas have been considered in detail. In the first model, the 
gas is distributed isothermally at the virial temperature of the mass. In the isentropic model, the adiabatic collapse 
of pre-heated gas onto an isothermal potential is considered, with the gas content of the halo constrained by the 
physics of adiabatic Bondi accretion. In both models we incorporate the Lacey & Cole (1993; 1994) distribution of 
halo formation times, and also the contribution of recombination radiation to the gas temperature and luminosity. 
The important conclusions are the following: 
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• Neglecting the distribution of halo formation times leads to a 10-20 per cent overestimate of the normalisation 
parameter as 

• Neglecting the contribution from recombination radiation has a significant effect on the slope of the luminosity- 
temperature relation, for kT < 4 keV. 

• In the isentropic model, we assume the entropy of the pre-infall gas has been increased by some, unspecified 
mechanism. For haloes with X-ray gas temperatures kT < 1 keV, the luminosity-temperature relation in this model 
is L oc T 5 . Both the slope and dispersion (due to the distribution of formation times) of this model are a good match 
to the available X-ray observations. The fit is improved, especially for fi = 1, if we assume a baryon fraction of 

n b /n a = o.i. 

• The isentropic model predicts that the gas temperature in haloes with M < 10 12 Mq will be significantly higher 
than the virial temperature; thus, making the assumption that the gas is at T v j r will lead to an overestimation of 
the mass by more than an order of magnitude. This discrepancy can be reduced somewhat if the gas accretion rate 
is lower than the critical adiabatic Bondi rate we have assumed here; however, we still predict that the gas fraction 
will increase proportionally to the halo mass in this regime. 

• The amount of heating required at z « 1 in the isentropic model is equivalent to heating halo gas at the 
turnaround epoch to 3 x 10 6 K for Q = 1 and 2 x 10 6 K for o = 0.3. 

• The slope of the observed luminosity function is well matched by the isentropic model for L < 10 43 ergs s _1 

• The isentropic model fails for more massive clusters, where the gas entropy becomes increasingly dominated by 
shocks. 

As better X-ray observations become available for less luminous, cooler groups of galaxies, the observed slope 
of the L-T relation in this regime will provide a strong test of the validity of the isentropic model. For the majority 
of clusters, the gas is clearly not isentropic or isothermal. These objects may be best described either by largely 



isothermal gas with a high entropy core ( Evrard, Henry 1991 ) , by the shock-heating model of Cavaliere et al. (1997) 



or by a model which treats the g uniformly mixed polytrope with 7a 1.2. 
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APPENDIX A: 



In this appendix we will discuss two different models in which pre-heated gas collapses adiabatically. We will make 
two assumptions about the behaviour of p g ,R and K a that are different from those which we make for our isentropic 
model, discussed in § 3.2.3. In both cases, we assume that p s ^b./pr = ^b/^o, as in the isothermal model. 

First, we can assume that K a is a universal constant, independent of halo mass. Then, from Equations 19 and 
|20| , the ratio of gas to dark matter mass within R vu will vary strongly with V c , in the sense that more massive 
clusters will have smaller gas fractions. The gas mass fraction in this case would vary by two orders of magnitude 
between haloes with V c — 300 km s _1 and those with V c = 1000 km s _1 . However, both observations and theory 
suggest that the baryon fraction in clusters may be greater than the universal mean, but do not deviate from this 
value by more than a factor of a few (White et al. 1992; White, Fabian 19*95] ; David, Jones, Forman 1995; Evrard 
1997). Thus, we reject this model as improbable on physical grounds. 

Alternatively, we can constrain the ratio of gas mass to dark matter mass to be a constant equal to the universal 
value, Mg/M = f2f,/f2 . Then from Equation |2C|, A=1.03 is independent of mass and Zf, which requires the entropy 
constant to scale as K Q oc V c 2 (from Equation 19). The entropy of the gas will still be constant for any given cluster, 
but will be higher for larger haloes; this requires a mechanism for pre-heating the gas that is more efficient for 
more massive haloes. In this case, since the parameter A is constant, T cm oc T g n ocV? oc T v ; r from Equations |2^ 
and |25| and, thus, the isothermal mass-temperature relation (T oc A/ 2 / 3 , Equation |l4|) still holds, apart from the 
normalisation. Furthermore, from Equation |26|, recalling that K oc V^, it holds that L oc V c 4 when bremsstrahlung 
radiation is the dominant emission mechanism and L oc when recombination radiation dominates; consequently, 
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L oc T and L oc T for the two cases, respectively. Interestingly, the above scaling relations are exactly the same as 



those derived in the isothermal case (§ 3.1). The relationships are not unique to the isothermal model. 
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